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ABSTRACT 



Context. Dwarf nova outbursts result from enhanced mass transport through the accretion disc of a cataclysmic variable system. 
Aims. We assess the question of whether these outbursts are caused by an enhanced mass transfer from the late-type main sequence 
star onto the white dwarf (so-called mass transfer instability model, MTI) or by a thermal instability in the accretion disc (disc insta- 
bility model, DIM). 

Methods. We compute non-LTE models and spectra of accretion discs in quiescence and outburst and construct spectral time se- 
quences for discs over a complete outburst cycle. We then compare our spectra to published optical spectroscopy of the dwarf nova 
SS Cygni. In particular, we investigate the hydrogen and helium line profiles that are turning from emission into absorption during the 
rise to outburst. 

Results. The evolution of the hydrogen and helium line profiles during the rise to outburst and decline clearly favour the disc- 
instability model. Our spectral model sequences allow us to distinguish inside-out and outside-in moving heating waves in the disc of 
SS Cygni, which can be related to symmetric and asymmetric outburst light curves, respectively. 

Key words. Accretion, accretion disks - Stars: dwarf novae - Novae, cataclysmic variables - Stars: individual: SS Cygni 



1. Introduction 

Dwarf novae (DN) belong to the non-magnetic cataclysmic vari- 
ables which are binary systems consisting of a white dwarf as 
primary component and an orbiting late-type main sequence star. 
Due to their close orbit, mass transfer from the secondary onto 
the primary via Roche lobe overflow occurs. Because of conser- 
vation of angular momentum an accretion disc forms around the 
white dwarf (Warner 1995). 

Dwarf novae are characterised by more or less regular out- 
bursts during which the system undergoes a rise in optical bright- 
ness of 2-6 magnitudes. The observed outbursts can be divided 
into two categories, depending whether the lightcurves are sym- 
metric or not. For asymmetric outbursts they are characterised 
by a fast rise and slower decline and show a delay of the rise 
in the UV against the optical. For symmetric outbursts UV and 
optical fluxes rise simultaneously on a longer timescale than the 
asymmetric outbursts. Sometimes for both types plateaus are ob- 
served during maximum. It is commonly accepted that the out- 
bursts are caused by a luminosity increase in the disc that arises 
from a temporarily increased mass transport through the disc. 
The origin of this increased mass transport, however, has been 
discussed controversially. 

According to the mass transfer instability model (MTI, Bath 
1975) an instability in the secondary star leads to a temporar- 
ily increased mass transfer rate M2 from the secondary onto the 
disc, so that also the mass transport through the disc increases. 
Because the instability originates in the secondary, the outbursts 
must start at the outer edge of the disc and proceed inwards in 
the framework of this model. 
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In the disc instability model (DIM, Osaki 1974) the mass 
transfer from the secondary is constant and the outbursts are at- 
tributed to thermal viscous instabilities in the disc that lead to a 
temporarily increased mass transport through the di sc. Meyer & 
Meyer-Hofmeister dl98ll) and Faulkner et al. d 19831) have shown 
that this instability is due to the local ionisation of hydrogen. 
Radial temperature and viscosity gradients lead to the propaga- 
tion of heating or cooling waves throughout the disc, which carry 
the whole disc over to outburst or quiescence, respectively. In 
particular this allows the outbursts to start any place in the disc 
so that outbursts can proceed inwards or outwards. 

Today the DIM is generally favoured over the MTI. This is 
due to the existence of a detailed theoretical framework for the 
DIM that can explain the outburst behaviour and the differen t 
outburst types in a natural manner. According to Smak (1984a) 
the different rise times for asymmetric and symmetric outbursts 
are caused by different propagation directions of the outbursts: 
for the asymmetric outbursts (type- A after Smak 1984a) the 
heating wave originates in the outer part of the disc and pro- 
ceeds inside moving with the mass stream, so that the heating 
wave can move relatively fast. Finally the hot inner part of the 
disc switches to outburst, so that the rise in the UV is delayed 
against the optical. In contrast to the symmetric type-B outbursts, 
the heating wave proceeds inside-out and therefore has to move 
against the mass stream, resulting in a relatively slow rise. In 
this case the hot inner parts of the disc are the first to switch to 
outburst so that optical and UV fluxes increase simultaneously. 

Furthermore, that no luminosity increase in the hotspot - the 
region of the disc, where the mass stream from the secondary 
impinges - is observed during an outburst contradicts the MTI 
model, because such a luminosity increase would be expected if 
the mass transfer from the secondary increases. 
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A possibility of distinguishing between both models by ob- 
servational data arises if one is able to decide whether the out- 
bursts proceed outside-in or inside-out. This can be achieved 
by comparing time-resolved spectra for an entire outburst cycle 
with the appropriate model spectra because quiescence and out- 
burst spectra differ significantly. In quiescence the optical spec- 
trum shows the strong hydrogen Balmer emission lines char- 
acteristic of an optically thin disc. In contrast, during outburst, 
broad absorption features in the Balmer series indicate an opti- 
cally thick disc. At the same time, the intensity in the blue wave- 
length range increases particularly strongly indicating a rise in 
the disc temperature. 

To this end we calculated time-resolved model spectra 
(Sect. 3) tailored to the dwarf nova SS Cyg. This is the bright- 
est known DN, showing an outburst brightness of V = 8.2 mag 
(Ritter & Kolb 2003), making it one of the best-studied DN. 
Before presenting these models in Sect. 3, we give a short 
overview of our approach in Sect. 2. The results are discussed 
in Sect. 4. 



2. Model assumptions 

To calculate the accretion disc spectra we use our accretion disc 
code AcDc (Nagel 2003; Nagel et al. 2004), which is based on 
the assumption of a geometrically thin disc (total disc thickness 
H is much smaller than the disc diameter). This allows us to 
decouple vertical and radial structures and, together with the as- 
sumption of axial symmetry, to separate the disc into concen- 
tric annuli of plane-parallel geometry. Then radiative transfer 
becomes a one-dimensional problem. 

Each of these disc rings, which are located at a given ra- 
dial distance r from the white dwarf, is assumed to be stationary. 
Thus it can be characterised by a constant mass transport rate M. 
The rate of energy generation from viscous shear then becomes 
independent of the kinematic viscosity Vk and can be parame- 
terised by the effective temperature 



T eft (r) 



3GM\M 



87nxr 3 
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1/4 



(1) 



(for example Warner 1995). Thereby Mi denotes the mass, r\ the 
radius of the primary white dwarf, G the gravitational constant, 
and cr the Stefan-Boltzmann constant. 

To get a self-consistent solution, the radiative transfer equa- 
tion, the hydrostatic and energy equilibrium equations, as well 
as the NLTE rate equations, that determine the occupation num- 
bers of the atomic levels, are solved simultaneously by an itera- 
tive scheme. Therefore detailed information about the involved 
atomic levels is needed, which is provided in the form of a 
model atom (cf. Rauch & Deetjen 2003). The kinematic viscos- 
ity, which is needed for the vertical structure calculation, c an b e 
parameterised by the a-approach of Shakura & Sunyaev d 1973b 



v k = ac s H 



(2) 



(where c s is the speed of sound and a < 1 a d imensionless pa- 
rameter) or after Lynden-Bell & Pringle ( 1 19741) by the Reynolds 
number Re 



a7 



Vgm 

Re 



(3) 



where is the Kepler velocity. We choose the latter approach 
that is numerically easier to implement because we save a further 
iteration to solve consistently for c s and H. 



Irradiation of the disc by the primary is considered via the 
upper boundary condition for the radiative transfer equation. For 
that purpose the irradiation angle /3 for each disc ring and the 
spectrum of the primary must be specified. The spectrum of the 
primary is parameterised by a blackbody temperature 7bb, or de- 
tailed white dwarf model atmosphere spectra are calculated. 

The complete set of input parameters, which we must pro- 
vide for each disc ring, thus consists of M\,ri,M,r,Re,f3,T\,t,. 
The spectrum of the complete disc is then obtained by integrat- 
ing the spectra of these disc rings for different inclination an- 
gles, the spectral lines are Doppler shifted according to the radial 
component of the Keplerian rotation velocity in the disc. 

As the accretion discs of dwarf novae are fed by a late-type 
main sequence star, we assume a disc composition of hydro- 
gen and helium with relative solar abundances. The model atoms 
used for the disc model calculations presented here contain the 
ionisation stages H i, H n, and He i-He m. The number of NLTE 
levels and lines considered is 15 and 105 for Hi, 29 and 61 for 
He i, and 14 and 78 for Hen. We consider the H opacity and 
Rayleigh scattering for H and He, which is important for the 
coolest regions of the disc model. In addition, the Ly„ line in 
the cool models for the quiescent disk is so broad that it con- 
tributes considerably to the source function in the optical band. 
The reason is that most of the hydrogen (about 90 - 99%) is neu- 
tral throughout most of the line-forming region (Fig. Q] second 
panel from top). 



3. Models 

In the following we present detailed models for the accretion 
disc of SS Cyg in outburst and quiescence. SS Cyg is the bright- 
est known DN and belongs to the U Gem type of DN. For 
our models we have chosen the orbital parameters according to 
Ritter & Kolb d2003l) . who give Mi = (1.19 ± 0.02) M for the 
mass of the white dwarf. According to the mass-radius relation, 
this corresponds to a white dwarf radius of 3.9 • 10 8 cm. Together 
with the mass of the companion M 2 - (0.704±0.002) M and the 
orbital period P = 6.6031 h, the tidal radius - the radius where 
the disc is disrupted by tidal interactions with the secondary - 
follows from 



Aidal = 0.60 ■ 



1 + q 



(4) 



(Hellier 2001) and amounts to r t id a i = 5.78 • 10 cm. Here a de- 
notes the distance between primary and secondary, and it can be 
calculated from the third Kepler law. q is the mass ratio j^. The 
minimal extension of the disc is given by the so-called circulari- 
sation radius 



1 + q 



(5) 



at which the angular momentum is equal to the angular mo- 
mentum at the Lagrange point Li. The Roche lobe and there- 
fore the distance of the Lagrange point L i mus t be calcu- 
lated numerically. After Plavec & Kratochvil (1964), however, 
for 0.1 < q < 10 the approximation 



r Ll = a- (0.500 -0.227 log q) 



(6) 



is possible. This finally leads to r c i rc = 1.65 • 10 10 cm. In this 
radial range we have increased the disc's outer edge r until 
the double-peaked line profiles matched the observation, so we 
chose r — 4 ■ 10 10 cm. A lower value would give line profiles 
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that are too broad due to the higher Kepler rotation velocity for 
smaller radii. 

The inner edge of the disc model was fixed from the fol- 
lowing arguments. Our model cannot be applied to the boundary 
layer expected at the transition from the disc to the primary, so 
we must truncate the disc well before the white dwarf. Despite 
the higher temperatures in the inner disc, this has little influence 
on the optical spectrum because the surface area of the inner 
parts of the disc is much smaller than that of the outer parts. Thus 
the inner parts can be neglected for modelling optical spectra. In 
contrast, in the UV range there will be strong imprints from the 
inner disc portions (cf. Fig.|2]i. 

For the inclination angle i we have chosen 40°, which is con- 
sistent with the value of i = (37 + 5) ° given by Ritter & Kolb 
(2003). All disc rings have been irradiated with a 50 000 K black- 
body spectrum. This temperature is compatible with the obser- 
vational results for the WD in SS Cyg (Long et al. 2005; Smak 
1984b). Tests with a 50 000 K white dwarf model atmosphere 
have shown that the blackbody approximation has little influence 
on the emerging disc spectra. The irradiation angle was set to 1°. 
According to the system geometry, this is possible but probably 
marks an upper limit. For such small angles, the irradiation in- 
creases the effective temperature of the disc ring compared to the 
value expected according to Eq.|2]only marginally. The relative 
differences are below 10~ 3 and decrease with increasing r. 

3.1. Outburst 

For the hot disc during outburst, we assume a constant mass 
transport rate through the disc of M = 4 • 10~ 9 M Q /yr and a 
constant viscosity of a x 0.30 according to the DIM. With those 
parameters, we calculated a disc model from 1 • 10 9 cm < r < 
4 ■ 10 10 cm by dividing the disc into 20 rings to obtain a smooth 
distribution of T^{r) with a maximal difference of ~ 3500 K 



70 t i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i y 




-6 -5 -4 -3 -2-10 1 

log m [g/cm 2 ] 

Fig. 1. Vertical structure of the hot disc at a distance of 7.35 ■ 
10 9 cm from the white dwarf. The physical variables are plotted 
against the column mass measured from the surface towards the 
midplane. 



Table 1. Parameters of the rings for the hot disc in SS Cyg. 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 



10 9 cm 



Re T cff [K] 



h 10 s cm 



1.00 
1.10 
1.22 
1.35 
1.50 
1.66 
1.84 
2.05 
2.30 
2.60 
2.97 
3.43 
4.02 
4.80 
5.85 
7.35 
9.65 
13.50 
21.00 
40.00 



3200 
3200 
3000 
2900 
2800 
2600 
2500 
2400 
2350 
2300 
2300 
2300 
2200 
2100 
1700 
1620 
1550 
1450 
1200 
500 



74912 
71056 
66935 
63013 
59074 
56742 
51916 
48403 
44873 
41350 
37798 
34259 
30705 
27135 
23610 
20080 
16525 
12969 
9404 
5862 



104 
110 
114 
119 
124 
125 
127 
128 
130 
131 
134 
139 
145 
155 
163 
192 
221 
258 
169 
1.61 



0.21 
0.23 
0.26 
0.30 
0.34 
0.39 
0.44 
0.50 
0.58 
0.67 
0.79 
0.93 
1.13 
1.40 
1.81 
2.38 
3.29 
4.93 
8.14 
16.18 



r to t is the total Rosseland optical depth from top to disc mid- 
plane and h = H/2 the vertical extension of the disc from the 
midplane. The other symbols are defined in the text. 



between neighbouring rings (see Table [TJ. The resulting disc is 
optically thick except for the outermost ring. This ring's effective 
temperature of 5862 K is more typical for a cold disc. 

As an example of the vertical structure of the disc model, the 
temperature, hydrogen ionisation fraction, density, and optical 
depth are plotted against the height above the disc midplane in 
Fig. Q] at a distance of 7.35 • 10 9 cm from the white dwarf. The 
temperature shows an inversion at the disc surface due to the 
heating by irradiation of the WD before a strong drop down at 
log m ~ -4 occurs and then the temperature rises slowly towards 
the disc's midplane. At log to * -1 the disc becomes optically 
thick. 

Figure [2] shows the integrated spectrum of the hot disc and 
the contribution of selected disc rings. Doppler broadening due 
to the Keplerian velocity is taken into account. In the optical it is 
characterised by hydrogen Balmer absorption lines. In the UV, 
strong absorption lines of the hydrogen Lyman series and He n 
appear. The latter lines originate in the inner disc rings, where 
T e g becomes high enough to populate He n levels. This high T e ff 
is also the reason the inner disc rings dominate the spectrum in 
the UV range despite their small surface area compared to the 
outer rings. The situation is completely different in the optical. 
There the disc spectrum is dominated by the outermost ring due 
to its large surface area, so the weak absorption line of He n at 
4686 A, which is visible in the inner ring spectra, is outshone by 
the much larger continuum flux of the outer rings. 

For a comparison to observational results, Fig. [3] shows our 
synthetic spectrum of the hot accretion disc with the spectra of 
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1000 2000 3000 4000 5000 6000 7000 



Fig. 2. Model spectra for the accretion disc of SS Cyg in outburst (uppermost curve). The other curves show the contribution of 
selected individual disc rings starting with the outermost ring (top) and then continuing towards the inner disc edge. The inclination 
angle is 40°. 
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Fig. 3. Observed spectra of SS Cyg (from Clarke et al. 1984) dur- 
ing rise to outburst (outburst maximum corresponds to 7/28/83). 
Overplotted is our synthetic spectrum for the accretion disc in 
outburst. The model flux was multiplied by a constant factor to 
match the observed continuum flux. 



Clarke et al. 0984), who observed SS Cyg during a rise to out- 
burst. 



3.2. Quiescence 



Wood et al. ([1986) studied the radial temperature distribution 
in the accretion disc of the DN Z Cha by eclipse mapping. In 
contrast to the T e ff oc r~ 3 ^ 4 power law expected for stationary 
accretion discs, they found a more or less constant value of the 
effective temperature at a level of several thousand Kelvin. This 
has been interpreted as a hint that the accretion discs of DN in 
quiescence are not stationary. Therefore we assumed a constant 



effective temperature of ~ 4200 K throughout the disc, which is 
in the typical range for cold discs. To achieve this temperature 
for all rings, we adjusted the mass transport rates through the 
rings (see Table |5J. We also adjusted the Reynolds number to 
get typical values for a in a disc in quiescence according to the 
DIM. The resulting kinematic viscosity is smaller than in the hot 
disc by a factor of 10, except for the outermost ring, where the 
viscosity is as high as for disc rings in outburst. For this disc 
ring, it was not possible to construct a low-viscosity model with 
strong emission lines. 



Table 2. Parameters for the ring models of the cold disc in SS 

Cyg. 



# 



lO'cm 



M[M /yr] 



Re 



h 



10 s cm 



4.00 
6.00 
8.00 
9.00 
10.00 
20.00 
40.00 



1.4 
4.0 
1.0 
1.4 
1.9 
1.4 
1.0 



10 12 

io- 12 

10-" 

io-" 
io-" 

10 -io 
10" 9 



19000 
16000 
13000 
13000 
10000 
3000 
500 



0.27 
0.24 
0.28 
0.30 
0.29 
0.28 
0.22 



0.71 
1.13 
1.59 
2.00 
2.34 
6.00 
14.56 



The number of disc rings required is much smaller than for 
the hot disc, as the change in spectral properties across the radius 
is marginal due to the constant temperature. Furthermore we did 
not extend the cold disc as far in as the hot disc, but truncated the 
model at r — 4 ■ 10 9 cm. This is again justified by the constant 
effective temperature, which prevents a strong contribution of 
the inner rings to the UV-flux in contrast to the case of the hot 
disc (see Fig.|4]l. 

The resulting disc spectrum (Fig. |4]i is com pared to an ob- 
served spectrum taken from Martinez-Pais et al. ( 1994) in Fig. [5] 
For that purpose the model spectrum has been normalised to the 
local continuum flux. In principle, the model reproduces the hy- 
drogen Balmer emission lines, but they are not as strong as in the 
observation. The He i emission lines are not seen in our model 
spectrum, and the reason may be that they form in the hot spot, 
which we have not included in our models. 
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Fig. 4. Model spectra for the accretion disc of SS Cyg in quiescence (uppermost curve). The other curves show the contribution of 
selected individual disc rings starting with the outermost ring (top) and then continuing towards the inner disc edge. The inclination 
is 40°. The spectral lines are getting broader due to the higher Kepler velocities of the inner disc rings. 
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Fig. 5. Observed spectra of the accretion disc of SS Cyg during 
quiescence (Martinez-Pais et al. 1994). Overplotted is our syn- 
thetic spectrum for the accretion disc of SS Cyg in quiescence 
(grey). The flux was normalised to the local continuum flux. 



As an example, the vertical structure of the cold disc is 
shown in Fig.|6]at a radial distance of 4.0 • 10 10 cm. In contrast to 
the hot disc, the temperature does not increase towards the disc 
midplane but declines monotonically. Towards the disc surface, 
the irradiation of the WD again causes a temperature inversion. 
The entire cold disc is optically thin. 



3.3. Rise to outburst 
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Fig. 6. Vertical structure of the cold disc at a distance of 4.0 • 
10 10 cm from the white dwarf. The physical variables are plotted 
against the column mass measured from the surface towards the 
midplane. 



To examine the spectral evolution from quiescence to outburst, 
we combined rings of the cold and hot discs to a sequence of 
disc models such that this sequence simulates the propagation of 
the heating wave throughout the disc. In this way we studied the 
two different cases of outside-in and inside-out moving heating 
waves to achieve further insight into the processes taking place 
in the disc during rise to outburst. For the outside-in outburst 
this sequence consists of five disc models in which the cold rings 



have been replaced by the next-neighbouring hot rings from the 
outside. The assembly of these models is shown in Table[3] 

The left panel of Fig. |7]shows the spectral evolution for this 
sequence from a pure cold disc to full outburst, as well as the 
left panel of Fig. |8]where the spectra are normalised to the local 
continuum. The Balmer series turns from emission to absorption 
immediately after the outermost disk ring flipped into the hot 
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state, because the overall disc flux in the optical is dominated by 
the flux of the outer rings and these are dominated by absorption 
in the hot state. 

Table 3. Assembly of disc models for the simulated outside-in 
outburst. 



# 


Step 1 


Step 2 


Step 3 


Step 4 


Step 5 


1 


4.00cold 


4.00cold 


4.00cold 


4.00cold 


4.00cold 


2 


6.00cold 


6.00cold 


6.00cold 


6.00cold 


6.00cold 


3 


8.00cold 


8.00cold 


8.00cold 


8.00cold 


7.35hot 


4 


9.00cold 


9.00cold 


9.00cold 


9.65hot 


9.65hot 


5 


lO.Ocold 


lO.Ocold 


13.5hot 


13.5hot 


13.5hot 


6 


20.0cold 


21.0hot 


21.0hot 


21.0hot 


21.0hot 


7 


40.0hot 


40.0hot 


40.0hot 


40.0hot 


40.0hot 



Numbers in "step" columns denote the radial position of the 
model in 10 9 cm, the following "cold" or "hot" label whether a 
cold or hot ring was used. 



Similarly we modelled the inside-out outburst by a sequence 
of five disc models. In the first step the cold disc's innermost 
ring is replaced by the hot rings that lie inside of its radial posi- 
tion. At the same time the disc is extended inwards to the inner 
boundary of the hot disc at 1 ■ 10 9 cm. For optical spectra, to 
which we will restrict our discussion in the following, this sim- 
plification can be justified because the inner disc rings will only 
contribute to the UV due to their high effective temperature and 
because they only cover a small surface area. In the subsequent 
steps, the next-neighbouring rings from the inside are replaced 
by hot ones. The complete assembly of the discs for the inside- 
out model sequence is shown in Table [4] 

Table 4. Assembly of disc models for the simulated inside-out 
outburst. 



# Step 1 Step 2 Step 3 Step 4 Step 5 



1 


l.OOhot 


l.OOhot 


l.OOhot 


l.OOhot 


l.OOhot 


2 


l.lOhot 


l.lOhot 


l.lOhot 


l.lOhot 


l.lOhot 


3 


1.22hot 


1.22hot 


1.22hot 


1.22hot 


1.22hot 


4 


1.35hot 


1.35hot 


1.35hot 


1.35hot 


1.35hot 


5 


1.50hot 


1.50hot 


1.50hot 


1.50hot 


1.50hot 


6 


1.66hot 


1.66hot 


1.66hot 


1.66hot 


1.66hot 


7 


1.84hot 


1.84hot 


1.84hot 


1.84hot 


1.84hot 


8 


2.05hot 


2.05hot 


2.05hot 


2.05hot 


2.05hot 


9 


2.30hot 


2.30hot 


2.30hot 


2.30hot 


2.30hot 


10 


2.60hot 


2.60hot 


2.60hot 


2.60hot 


2.60hot 


11 


2.97hot 


2.97hot 


2.97hot 


2.97hot 


2.97hot 


12 


3.43hot 


3.43hot 


3.43hot 


3.43hot 


3.43hot 


13 


4.02hot 


4.02hot 


4.02hot 


4.02hot 


4.02hot 


14 


4.80hot 


4.80hot 


4.80hot 


4.80hot 


4.80hot 


15 


5.85hot 


5.85hot 


5.85hot 


5.85hot 


5.85hot 


16 


6.00cold 


7.35hot 


7.35hot 


7.35hot 


7.35hot 


17 


8.00cold 


8.00cold 


9.00cold 


9.65hot 


9.65hot 


18 


9.00cold 


9.00cold 


lO.Ocold 


lO.Ocold 


13.5hot 


19 


lO.Ocold 


lO.Ocold 


20.0cold 


20.0cold 


21.0hot 


20 


20.0cold 


20.0cold 


40.0cold 


40.0cold 


40.0cold 


21 


40.0cold 


40.0cold 









The right panels of Figs.|7]and[8]show the spectral evolution 
for this sequence from a pure cold disc to full outburst. In con- 



trast to the outside-in scenario the hydrogen Balmer emission 
only diminishes slowly during rise to outburst, while increasing 
absorption wings appear. 

One has to keep in mind that the steps of our sequences are 
not equidistant in time. For an outside-in outburst, for exam- 
ple, the heating wave moves inwards quickly; and according to 
our models, the spectral lines change from emission to absorp- 
tion as soon as a part of the outer region is in outburst, so one 
will observe an absorption-line spectrum during most of the rise 
of an outside-in outburst. In the case of an inside-out outburst, 
the heating wave moves rather slowly outwards. As our models 
show, the line spectrum of the disk does not change to pure ab- 
sorption until the outermost regions are in outburst. Hence one 
would observe an emission-line spectrum most of the time of an 
inside-out outburst. 

3.4. Decline 

Decline from outburst to quiescence is mediated by a cooling 
wave. T his wa ve always propagates outside-in as argued by e.g. 
Warner ( 1995), so that the cooling can be studied by examining 
the inside-out moving heating wave of Table [4] in reverse order. 
According to the right panel of Fig. [7] the hydrogen Balmer lines 
evolve smoothly from absorption to pure emission. 

4. Results and discussion 

We used the models presented above to determine the nature of 
the outbursts in SS Cyg by comparing them to spectra available 
in the literature. This turned out to be rather difficult, because 
adequate time-resolved spectra during rise to outburst are very 
rare even for the well-studied case of SS Cyg. 

Martinez-Pais et al. d 19961) presented time-resolved spectra 
for different outbursts of SS Cyg. Among them are two spec- 
tra taken during a rise to an outburst of the symmetric type. 
According to their Fig. 2 the hydrogen Balmer emission lines 
decrease slowly between these two spectra, while the absorp- 
tion wings increase like in the right panel of our Fig. [7] for the 
inside-out outburst. This leads to an identification of symmetric 
outbursts in SS Cyg with inside-out outbursts, which is in good 
agreement to Smak's (1984a) description of type-B outbursts, 
but contrasts to the original conclusion of Martinez-Pais et al. 
d!996l) . They interpret the late appearance of the Hen 4686 A 
line as a consequence of an outside-in propagating heating wave, 
because they assumed that the Hen 4686 A line originates in 
the hot inner part of the disc. This is questionable in light of 
our models. The line should be significantly broader due to the 
higher Kepler velocity if this assumption is true. 

"Some indication of an increase in the hot spot's vertical size 
and, perhaps, brightness" leads them to conclude further that 
symmetric outbursts are connected with an instability of the sec- 
ondary star, so they favoured the MTI for this outburst. This is 
put into question by our models, which indicate an inside-out 
outburst in favour of the DIM. 

The observations of Clarke et al. (1984|) cover a complete 
outburst cycle of SS Cyg. The outburst shows an asymmetric 
light curve. Three of its spectra before and during rise, as well 
as during maximum, are shown in Fig. [3] The spectral evolu- 
tion of this outburst differs significantly from what was observed 
by Martinez-Pais et al. d!996l) : quiescent Balmer line emission 
abruptly disappears at the onset of rise to maximum, before full 
absorption sets in during the rise. This fits our outside-in model 
sequence, indicating that the asymmetric outbursts are indeed 
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Fig. 7. Spectral evolution between 4000 and 5300 A for an outside-in outburst (left panel) and an inside-out outburst (right panel). 
The lowermost graphs show a pure cold disc that evolves to full outburst (uppermost graph) over steps 1 to 5 of Tables [3] and H] 
respectively. 




4000 

MA] 



Fig. 8. Spectral evolution between 4000 and 5300 A for an outside-in outburst (left panel) and an inside-out outburst (right panel). 
The spectra are normalised to the local continuum. The lowermost graphs show a pure cold disc that evolves to full outburst 
(uppermost graph) over steps 1 to 5 of Tables[3]and|4j respectively. 



connec ted to outside-in, i.e. type-A outbursts following Smak 
(Il984al) . 

The fact that the observed spectra of Clarke et al. ( 1984) dur- 
ing full maximum do not show pure absorption lines like our 
model might be a consequence of the radial extension of the disc 
during outburst. This is not considered in our model, although 
it is expected in the DIM due to the higher transport of angular 



momentum during outburst. It might be possible that this por- 
tion of the disc outshines the basic inner part of the disc due to 
its large surface area. If this portion of the disc then has com- 
parable properties to the current outermost grid point, which is 
relatively cool and emission-dominated, or at least continuum- 
dominated, the resulting model spectrum of the extended disc 
might show no absorption lines anymore. 
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For decline from outburst to quiescence, another study by 
Hessman et al. ( 1984|) exists. Their Fig. 5 is in good agreement 
with the right panel of our Fig. [7] if read from top to bottom, 
which means that they witnessed an outside-in propagating cool- 
ing wave. This again agrees with the DIM. 

During decline the 4686 A line of He n often shows a promi- 
nent emission feature like in the study of Hessman et al. 41984). 
This does not appear in our model. If it originates in the disc, it 
must arise from the inner parts, because only there does T e ff be- 
come high enough to populate Hen levels. Accordingly, the in- 
ner rings of our hot disc model show He n 4686 A, however, not 
in emission but in absorption. This supports the results of Unda- 
Sanzana et al. (2006), which identified the gas stream/disc im- 
pact region as the origin of the He n 4686 A emission by means 
of Doppler tomography. 

That Hessman et al. (1984) observed central emission peaks 
in the Balmer lines right from the beginning of the decline might 
be indicating that not the complete disc but only the inner parts 
participated in the outburst. If for example the disc stays cold for 
r > 10 • 10 10 cm, we only have to compare the lower five curves 
of t he rig ht panel in Fig. [7] to the observation. As Hessman et 
al. l[l984) observed a short outburst without plateau, this would 
be again in good agreement with the DIM. There the short out- 
bursts are attributed to discs that are not completely in outburst, 
while plateaus are supposed to appear if matter is accreted with 
a constant rate through a disc that is completely in outburst. 

It will be interesting to extend our study to the UV range 
where especially the C iv 1550 A line shows a similar behaviour 
to the Balmer lines. For that purpose, heavier elements must be 
included in the model calculations, and the influence of the disc 
wind, which becomes obvious in the P Cyg shaped profile of the 
C iv 1550 A line in outburst, will be considered. 

Another question is the influence of metal opacities, which 
we have neglected here, on the hydrogen and helium lines. From 
our experience in working on stellar atmospheres (O and sdO 
stars), we would predict that metal line blanketing and surface 
cooling will produce slightly deeper H and He absorption lines 
in the hot disk; however, we expect no qualitative change in the 
optical spectrum. The situation is different in the cool ring mod- 
els, which are optically thin. Work is in progress to investigate 
the metal line blanketing problem. 
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5. Summary 

In this paper, we have presented NLTE model calculations for 
the accretion disc of SS Cyg in outburst and quiescence. The 
resulting synthetic spectra describe the observed optical spec- 
tra and their transition well from absorption during outburst to 
emission during quiescence. 

Simulations of the spectral evolution for outside-in and 
inside-out propagating heating waves were carried out. We com- 
pared them to published observations and conclude that symmet- 
ric outbursts belong to the inside-out type. This confirms DIM 
expectations (e.g. Smak 1984a), which are based on rise-time 
arguments and explicitly excludes the MTI model. In contrast, 
asymmetric outbursts seem to be outside-in outbursts. 
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